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We investigate the cosmological efTects of a neutrino interaction with cold dark matter. We 
postulate a neutrino that interacts with a "neutrino interacting dark matter" (NIDM) particle 
with an elastic-scattering cross section that either decreases with temperature as or remains 
constant with temperature. The neutrino-dark-matter interaction results in a neutrino-dark-matter 
fluid with pressure, and this pressure results in diffusion-damped oscillations in the matter power 
spectrum, analogous to the acoustic oscillations in the baryon-photon fluid. We discuss the bounds 
from the Sloan Digital Sky Survey on the NIDM opacity (ratio of cross section to NIDM-particle 
mass) and compare with the constraint from observation of neutrinos from supernova 1987A. If 
only a fraction of the dark matter interacts with neutrinos, then NIDM oscillations may affect 
current cosmological constraints from measurements of galaxy clustering. We discuss how detection 
of NIDM oscillations would suggest a particle-antiparticle asymmetry in the dark-matter sector. 
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I. INTRODUCTION 

Flat cosmological models with baryons (about 5% 
of the total energy content of the universe) , cold dark 
matter (CDM, 25%), cosmological constant (or dark 
energy, 70%) and an adiabatic, nearly scale-invariant 
spectrum of density fluctuations explain most cosmo- 
logical observations. However, we still lack a satis- 
factory understanding of both dark matter and dark 
energy, a puzzle for both particle physics and cosmol- 
ogy- 

The most favored candidates for dark matter are 
cold, coUisionless massive particles, which are non- 
relativistic for most of the history of the Universe 
and so can cluster gravitationally during matter dom- 
ination. Candidates for these dark-matter parti- 
cles can be found in supersymmetric extensions of 
the standard electroweak model — namely neutralinos 
with mass on the order of 100 GeV Q — or in other 
theories (e.g., the axion, which may arise in the Peccei- 
Quinn mechanism 0). These cold-dark- matter mod- 
els account well for cosmic microwave background 
(CMB) observations on the largest scales as well as 
measurements of the large-scale distribution of galax- 
ies. 

However, observations on galactic and sub-galactic 
scales may conflict with the predictions, from numer- 
ical simulations and analytic calculations, of CDM 
models. Indeed, cold and coUisionless dark-matter 
models seem to predict an excess of small-scale struc- 
tures 0, and numerical simulations 3 predict far 
more satellite galaxies in the Milky Way halo than 
are observed. 

Several solutions have been proposed to explain 
these discrepancies, for exaniple, inflationary models 
with broken scale invariance However, most other 
explanations invoke modifications to the properties of 
dark-matter particles. For example, a warm-dark- 
matter candidate, like a sterile neutrino, has been 



suggested because it suffers free streaming and sup- 
presses the matter power spectrum on small scales Q • 
A dark-matter particle that results from decay of a 
short-lived charged particle can also suppress small- 
scale power i7|. Other possibilities include a dark- 
matter particle that interacts with others particles 
such as photons H, 0, or neutrinos 0, 0| or self- 
interacting dark matter For ^ review of differ- 
ent alternative scenarios to standard coUisionless cold 
dark matter, see Ref. [T^ . 

In this paper, we investigate the possibility of a 
neutrino interacting dark matter (NIDM) component. 
If dark matter and neutrinos interact, there was an 
epoch in the very early universe in which they were 
strongly coupled. Dark-matter perturbations that en- 
tered the horizon during this period would then be 
erased because of diffusion damping, and the suppres- 
sion scale will depend on the dark-matter-neutrino 
interaction. Even if only a fraction of the dark mat- 
ter interacts with neutrinos, a pattern of oscillations 
in the matter power spectrum arises, much like the 
oscillations in the baryon-photon fluid. 

In the following, we limit our study of DM-neutrino 
couplings to effects on cosmological scales in the fre- 
quency range smaller than k < 0.2/iMpc~^ — i.e., on 
scales where linear perturbation theory is viable. We 
consider flat cosmological models with an adiabatic 
and nearly scale invariant spectrum P{k) ~ k" of 
density perturbations where n = 0.97. Unless explic- 
itly stated, the energy content of the Universe corre- 
sponds to the standard ACDM model with baryons 
contributing as Qb — 0.05 and a cold-dark-matter en- 
ergy density Qdm — 0.25. We also choose a Hubble 
parameter h = 0.73, a standard value 3.04 for the ef- 
fective number of (massless) neutrinos ^^'^ dark 
energy in the form of a cosmological constant with 
equation-of-state parameter w = —1. The interaction 
between dark matter and neutrinos is given in terms 
of the opacity Q ~ (cTdm-i/I^^D/^dm, the ratio of the 
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thermal averaged dark-matter-neutrino cross section 
to the mass of the dark-matter particle. 

The paper is organized as follows. In the next Sec- 
tion, we discuss a class of models of neutrino-dark- 
matter interaction for both scalar and spinor dark- 
matter candidates and obtain an estimate for the 
opacity Q. In Section III, we outline the cosmological 
consequences of a NIDM component, and we compare 
those predictions with the latest data on galaxy clus- 
tering from the Sloan Digital Sky Survey (SDSS). In 
Section IV, we consider astrophysical constraints, par- 
ticularly those from observation of neutrinos from su- 
pernova 1987A. Finally, in the last Section, we report 
our conclusions. 



II. THE NEUTRINO-DARK-MATTER 
INTERACTION 

The possibility of new neutrino interactions with ex- 
otic matter fields and their cosmological implications 
have been recently considered by many authors. Cou- 
plings with a light scalar or pseudoscalar boson, as in 
the Majoron model [H IH iH 13 , can have sizeable 
effects on the CMB and the power spectrum of large- 
scale structure (LSS), and might lead to a neutrinoless 
Universe for a massless or very light scalar field 19] . 

A different class of interactions of neutrinos with 
DM, and more generally of DM with the electro- 
magnetic plasma as well, was put forward in a se- 
ries of papers discussing the possibility that DM con- 
sists of particles with mass TOdm in the MeV range 
I2D, HI In in Hi- If the rehc abundance of these 
particle is produced via the freezeout of annihilation 
processes and yet corresponds to the observed DM en- 
ergy density today, they must interact with stronger 
than weak interactions, since the cross section for an- 
nihilations via the exchange of a massive particle as, 
e.g., a vector boson, typically decreases at the freeze- 
out as the square of the DM-particle mass. 

The implications of this scenario for big bang nucle- 
osynthesis (BBN) have been considered in Ref. p^ . 
With DM-neutrino interactions producing the re- 
quired DM annihilation cross section of a few pico- 
barns (for an s-wave annihilation) at decoupling, the 
mass range mdm < 10 MeV is disfavored. In fact, 
tightly-coupled DM particles give in this case a non- 
negligible contribution to the total energy density, and 
the neutrino-photon temperature ratio is increased be- 
cause of the entropy release from DM to neutrinos. 
Both effects conspire to produce an order one extra 
effective neutrino species and thus an increased pri- 
mordial helium abundance, which is difficult to recon- 
cile with the present determination of the ^He mass 
fraction; see, e.g., Refs. jM^?, 28, 29]. 

Neutrino coupling to DM particles with mass in 
the MeV range can also be bound using neutrino 
fluxes from type II supernovae [s^l, such as SN1987A, 
again resulting in a constraint mdm > 10 MeV. It is 



worth observing however, that this bound only ap- 
plies if these particles have large elastic scattering 
cross sections with nucleons, larger than neutrino- 
nucleon ones, whose effect is to shift outward the neu- 
trinosphere and thus leads to a lower neutrino decou- 
pling temperature. 

In the following, we consider a DM-neutrino cou- 
pling even stronger than those considered in Refs. [20I 
I21I m m 1^ , so that the resulting large scattering 
cross sections might lead to detectable effects in the 
LSS power spectrum. In fact, as already emphasized 
in Ref. [23|, a DM-neutrino coupling which leads to 
an order pico-barn annihilation cross section can only 
affect the power spectrum on very small scales, since 
scattering processes freeze out quite early in time. In 
this case, LSS forms as in the presence of the usual 
collisionless DM fluid, while neutrinos act via free 
streaming and suppress power on scales smaller than 
the acoustic horizon at the time they become non- 
relativistic. 

We assume mdm > 10 MeV, since our large DM- 
neutrino coupling implies that the bound obtained in 
Ref. [2^ applies a fortiori Notice that our working 
hypothesis potentially leads to the problem of how 
a sizeable relic abundance of DM particle actually 
forms, since a stronger coupling with neutrinos en- 
hances their annihilation rate and leads to negligi- 
ble relic abundance today. We will comment on this 
point later and emphasize that our scenario requires 
an asymmetry between DM particle and antiparticle 
produced at an early stage. 

The interaction Lagrangian density depends upon 
the spin content of the DM field, denoted in the fol- 
lowing by ip. We consider here the case of a non- 
self-conjugated scalar particle or a Dirac spinor and 
interaction terms that admit a conserved global U{1) 
charge with transforming with a phase under the 
corresponding transformations. For a scalar ip we can 
write 

Ant = hFRVLi^ + h.c. , (1) 

where F is a spinor field. Similarly, DM coupling to 
neutrinos can be also introduced via the interaction 
with an intermediate vector-boson field J7^, 

+ 9vVl1^vlU^ . (2) 

In both cases, we consider the F 01 U field to have 
masses of order of MeV or larger to forbid -0 decays 
at tree level, which would erase any relic abundance 
of V' unless the couplings are tuned to be very small. 
In a way, this assumption is far from being ad hoc for 
the case of a coupling as in Eq. . If the F field is 
lighter than ■0, we simply have to shift our perspective 
and consider F as a candidate for light dark matter 
rather than the ■0 held. 

In the range of neutrino temperature T < MeV we 
are interested in, the thermally averaged ^-neutrino 
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scattering cross section is 



{ml - m2^)2 ' 



(3) 



for an F fermion exchange, unless the and F fields 
are degenerate in mass, in which case the low-energy 
transfer scattering cross section gets the usual Thom- 
son behavior, 



^dm — u 



\h\' 



dm 



Similarly, for a U coupling we have 
{cTdm^i'lvl) ^ glgl — 



(4) 



(5) 



The Lagrangian density, Eq. Q), also describes the in- 
teraction of (chiral) fermionic dark matter with neu- 
trinos via the exchange of a massive scalar field, with 
the obvious redefinition ■0 <-> f , and scattering cross 
sections are again given in this case by Eqs. © or Q). 

Finally, interaction of a Dirac DM field with neu- 
trinos via a vector-boson interaction Lagrangian, 

+ 9i,VLl^VLUt, , (6) 



gives, for rudm > T, 

{(Jdm-v\v\) ^ g^gl (c| + c 



rp2 

clCr) — . (7) 



A light U boson with an order MeV mass coupled 
to charged leptons might affect the electron-neutrino 
scattering cross section at low energy, while mea- 
surements shows no significant deviations from the 
standard electroweak-model result [sil l3^ (see also 
Ref. 21] for a detailed analysis on this issue). A pos- 
sible way out is of course to suppress the value of the U 
coupling to ordinary matter, including neutrinos, with 
respect to its coupling to DM particles. In this case, 
however, neutrino-dark- matter scattering would be 
quite small and again no observable effects on LSS 
can be obtained. Another possibility, though less ap- 
pealing, is to assume that the U boson couples mainly 
to neutrinos and very weakly to charged leptons. 

Regardless of the particular nature of DM parti- 
cles and the particular coupling to neutrinos, provided 
their mass as well as the mass of the exchanged par- 
ticle is in the range of MeV or larger, we see from our 
discussion that the typical thermally averaged scatter- 
ing cross section with neutrinos for T < MeV has two 
possible distinct behaviors, either decreasing as or 
constant for mass degeneracy of DM and intermediate 
scalar/fermion particle F. In particular, it is useful to 
define the DM-neutrino opacity, the thermally aver- 
aged scattering cross section over DM mass ratio, as 
follows 



{ad7n-^\v\) _ 



rridm 



(8) 



for the behavior as in Eqs. Q, ©, or 0. Here, 
a denotes the scale factor, normalized to unity at the 
present time. 

The opacity can be written in terms of the DM- 
neutrino coupling, which will be generically denoted 
by g and a mass scale M of order MeV or larger. For 
example, for a scalar DM particle coupled to neutrinos 
via the exchange of a fermion particle F we have g = 



\h\ and 



m 



dm 



[see Eq. while for a U 



exchange, g'^ = g^g y and — fnij [see Eq. (jSJ]. In 
this notation, and using the known value of neutrino 
temperature today, we get 



Q2 



g 



1 



(M/MeV)4 mdm/MeV 



■lO^^'^cm^ MeV"^ (9) 



Similarly, for the case of a constant scattering cross 
section as in Eq. we define 



(pdm-v\v\) _ 

= WO 



mdn 



which gives 



Qo 



{mdm/MeVy 



lO-^^cm^ MeV" 



(10) 



(11) 



Upper bounds on these parameters will be discussed 
in the following Section using the LSS power spectrum 
as well as other astrophysical constraints. 

We now come back to the issue of the relic abun- 
dance of DM particles ip. In the usual scenario, the 
present value of the DM energy density results from 
the freezeout of annihilation processes at tempera- 
tures of the order of mdm/20; see, e.g., Ref. [33 . As 
we will see in the next Section, the key parameter 
entering the Euler equation ruling the DM fluid per- 
turbation is the effective DM-neutrino scattering rate, 
defined by 



Pdrt 



~^dm^dm — iy : 



(12) 



with Udm the DM number density and and pdm 
the energy density of neutrinos and DM particles, re- 
spectively. If we assume that rridm > 10 MeV and 
that scattering processes with neutrinos are still ef- 
fective for small temperatures T ^MeV, thus leaving 
an imprint on the LSS power spectrum and CMB, an- 
nihilation to neutrinos would reduce the DM energy 
density to a very tiny value today. 

To show this let us assume that indeed the annihi- 
lation processes iptp vv freeze out at temperature 
of the order of Ty — mdm/^O^. Wc first consider the 



^ In general, DM particles are coupled to the electromagnetic 
plasma as well and annihilate into e+e" pairs. In this case, 
the corresponding annihilations will be assumed to freeze out 
at Tf or earlier. 
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case of a (non-self-conjugated) scalar DM field. For 
the coupling to neutrinos via the exchange of a mas- 
sive fermion F [see Eq. ([^l], we get for the annihilation 
cross section, 



'dm 



T 



(13) 



Notice that the s-wave annihilation contribution is 
suppressed by the square of neutrino mass, and has 
been neglected. This is a consequence of the fact 
that the coupHng is chiral; see, e.g., Ref. 0. Using 
Eq. (O, we see that at Tf the ratio of the thermally 
averaged neutrino scattering rate over the annihilation 
rate Tann = <yann\v\ Udm is of the Order of 



m 



dm J 



Tf 



^ann (ml - m^J' 

when mp ^ mdm [see Eq. Q] or 

Esc . I^dm 



rridn 



f 



(14) 



(15) 



when mp = mdm [see Eq. (0J]. This means that if 
annihilation freezes out at T ~ mdm-, neutrino scat- 
terings on DM particles are also largely ineffective at 
temperature smaller than MeV. In this case, no signa- 
tures of neutrino-DM interactions can be constrained 
by present observations. The largest scale where DM- 
neutrino interactions can leave an imprint corresponds 
to a wavenumber 



27r H{z,c) 

C 1 + Zsc 



• IQ-^hMpc-^ , (16) 



eq 



where H{zsc) is the Hubble parameter at the redshift 
Zsc where DM-i/ interactions freeze out and Zf,q is the 
redshift of matter-radiation equality. For Zsc ^ 10^, 
corresponding to T 0.01 MeV, we obtain k ^ lO'^h 
Mpc"^ 

Let us now consider the case of a C/ -exchange inter- 
action. The right order of magnitude (10 pb) of the {p- 
wave) annihilation cross section requires a light U bo- 
son mass, of the order of MeV, and coupling h ~ 10"'^. 
For T < mdm, we get 



'dm 



T 



^ '^'-^ [ml - AmlJ^ mdr, 



(17) 



Comparing this result with Eq. (jSj), we get at T/, 



{ml-AmlJ^ Tf 
mfr mdn 



(18) 



Again neutrino scatterings freeze out quite early 
around Tf. 

For spin-1/2 DM (Dirac or Majorana) particles and 
F or U couplings of Eqs. ^ and ©, one can reason 
along the same lines, with similar results. 



The case of a neutral scalar DM particle ip is quite 
different. In this case in fact, the F coupling to neutri- 
nos corresponds to an annihilation cross section which 
vanishes in the limit of massless neutrinos. On the 
other hand, couplings of DM to charged leptons can 
produce the correct relic abundance ii mp is of the 
order of 100 GeV to 1 TeV 21]. However, for such 
large values of mi?, the scattering cross section is very 
small too (in particular it vanishes in the local limit 
mp ^ oo due to a cancellation between the s- and 
M-channel amplitudes), so the case of self-conjugated 
scalar-DM-neutrino coupling is of no interest for the 
purposes of this paper, and is not expected to produce 
any observable features in the LSS. 

Summarizing, if we assume that DM couples to 
neutrinos strongly enough to produce observable ef- 
fects that can be constrained by CMB and LSS ob- 
servations, we have to abandon the idea that relic 
DM density formed via the usual mechanism based 
on freezing of DM annihilation processes at tempera- 
tures T mdm- 

How can our scenario be reconciled with the ob- 
served DM contribution to the present energy density 
of the Universe? First of all, it should be mentioned 
that only a fraction of the total DM could be cou- 
pled to neutrinos. In fact, we will consider this case 
too in the following. However, this would not repre- 
sent a solution to the problem, since this neutrino- 
coupled component would completely annihilate into 
neutrinos at Tf. The more interesting possibility is 
therefore that there is a particle-antiparticle asymme- 
try produced at higher temperatures in the DM sector 
coupled to neutrinos, very much like the mechanism 
by which the baryon (and lepton, in the framework 
of leptogenesis) number is produced in the early Uni- 
verse. Indeed, this possibility is also motivated by 
the intriguing observation that the parameters f2f, and 
rirfm only differ by a factor five today, yet their pro- 
duction mechanism is usually considered to be quite 
distinct, with very few exceptions 13 4. [351 : see also 
the discussion on this point in Ref. j^J. Though a 
more detailed theoretical analysis of this possibility is 
perhaps still needed, we think that this idea still rep- 
resents a stimulating possible scenario. In this case, 
it is meaningful to check to what extent present data 
can constrain strong couplings of DM particles with 
mass in the range MeV or higher to neutrinos. This 
scenario requires that DM particles are excitations of 
a non-self-conjugated field, such as a complex scalar 
field or a Dirac field, and that a particle-antiparticle 
asymmetry in the DM sector 77^ = (n^ — n:^) / n-y has 
been produced at some early stage, of the order of 



b mdm 



(19) 



with mp the proton mass and rjB the baryon-to-photon 
ratio, TjB ~ 6.3 x 10~^°; see e.g., Ref. [23. The analy- 
sis performed in the following Sections relies on these 
assumptions. 
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We conclude this Section by reporting the expres- 
sion of the smaller-wave mode k for which we expect to 
see the effects of DM-i^ scatterings in the dark-matter 
perturbation in terms of the opacities Q2 and Qq in- 
troduced in Eqs. Q and ifTTl) . This might be useful 
to understand the results reported in the following. 
From the definition of z^c. 



1e+06 



(20) 



and using standard values for the neutrino tempera- 
ture today and assuming no extra relativistic degrees 
of freedom in addition to photons and neutrinos, it is 
easy to get from Eq. II6II 

/c~0.2( j /iMpc~\ (21) 

or, for a constant scattering cross section, 
„ ,'10-22gm2 MeV 

fc - 0.2 X 10" 



/iMpc"^ 
(22) 



III. NIDM AND STRUCTURE FORMATION 

In order to include a neutrino-dark-matter interac- 
tion, we modify the standard Euler equations resulting 
in (using conformal time) , 



-—fdm 



+ andm(Jdm-v{Ov — ddm) , (23) 

oPdm 



dTldYn^dm — uiPdn 



(24) 



where momentum conservation in scattering processes 
has been accounted for. With 9 we denote the velocity 
perturbations, the subscripts "dm" and 'V" standing 
for neutrino-interacting dark matter and neutrinos, 
respectively. The quantity andmO'dm-w is the differen- 
tial opacity and gives the scattering rate of neutrinos 
by dark matter. 

As in the case of the baryon-photon interaction, we 
neglect both the shear term fc^cr and the term c^„fc^(5, 
where Cdm is the sound speed of the dark-matter fluid. 

As seen in the previous Section, the coupling be- 
tween neutrinos and dark matter can be parameter- 
ized through a cross section that either decreases as 
a~^, or takes a constant value. In these two cases, the 
parameters determining the DM perturbations are the 
opacities (cross section to DM-mass ratios) Q2 and Qq 
defined in Eqs. © and ifTUIl . 

First consider the case {(Jdm-v\''^\) oc a"^ (similar 
considerations can be made for the constant cross sec- 
tion). In Fig. ^ we show what happens when a per- 
turbation of wavenumber k = 1.04 /i Mpc~^ enters the 
horizon for different values of Q2. If the coupling is 




1e-05 



ie-04 0.001 0.01 
scale factor a 



FIG. 1: Dark-matter perturbations of = 1.04 /i Mpc~^; 
the opacity Q2 is in unit of cm'^ MeV~^. Damped oscilla- 
tions are clearly seen for Q2 = lO""*^ cm^ MeV~^. 



zero, we have the standard picture. The mode en- 
ters the horizon in the radiation-dominated era, and 
it starts to grow first logarithmically and then lin- 
early with the expansion factor (during matter dom- 
ination). When the same mode enters the horizon 
with Q2 = 5 X 10"'^"' cm^ MeV"\ the growth is 
nearly zero during the radiation epoch, while the mode 
starts growing linearly with the scale factor during 
matter domination, since the coupling with neutrinos 
becomes negligible in this stage for the chosen value 
of Q2. 

The situation is different when we consider a 
stronger coupling, say Q2 = 10~^^ cm^ MeV~^. When 
the perturbation enters the horizon, dark matter is 
coupled with neutrinos and this results in a series of 
oscillations until decoupling is reached. Notice that 
the amplitude of oscillations decreases near decou- 
pling, because the decoupling itself is not instanta- 
neous and so we see diffusion damping for the dark- 
matter-neutrino fluid. In Fig. 121 we plot several mat- 
ter power spectra for different values of the dark- 
matter-neutrino interaction for both the couplings 
considered. The effect of the dark-matter-neutrino 
interaction can be seen on small scales in the mat- 
ter power spectrum. Larger couplings will correspond 
to later epochs of neutrino-DM decoupling and to 
a damped oscillating regime on larger scales. Fi- 
nally, in Fig. 131 we plot the angular power spectra 
of CMB anisotropics for two models with and without 
DM-neutrino coupling. For the value of Q2 we con- 
sider, which is already at odds with current clustering 
data, there is a small enhancement in the small-scale 
CMB anisotropics. The reason for this is that the 
anisotropic stress in the neutrino relativistic compo- 
nent is reduced due to the coupling. In other words, 
neutrinos are no more a fluid with a "viscosity pa- 
rameter" cjjg = 1/3. This parameter, introduced 
in Ref. [s^, controls the relationship between veloc- 
ity/metric shear and anisotropic stresses in the neu- 
trino background [33 ■ The value of c^^^ will be close 
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FIG. 2: Several matter power spectra with different opac- 
ities Q2 (top panel) and Qo (bottom panel) between 
dark matter and neutrinos; Q2 and Qo are in units of 
cm^ MeV~\ 




10000 



FIG. 3: Angular power spectra with and without dark- 
matter-neutrino coupling. A small enhancement (~ 10%) 
of the height of the peaks on small scales is observed. 



to c^j^ — and this implies a small enhancement of 
the small-scale peaks, at the level of 10% In 
order to bound the strength of DM coupling to neu- 
trinos, we consider the real-space power spectrum of 
galaxies in the Sloan Digital Sky Survey (SDSS) us- 
ing the data and window functions of the analysis of 
Ref. 113. To compute the Hkelihood function C^^^s 



for the SDSS, we restrict the analysis to a range of 
scales over which the fluctuations are assumed to be in 
the linear regime {k < 0.2/i~^Mpc), and we marginal- 
ize over a bias b considered to be an additional free 
parameter. Since no relevant signature is expected 
on CMB anisotropics for the values of Q2 and Qo we 
consider, we do not include CMB-anisotropy data in 
the analysis and assume a cosmological concordance 
model with Cl\ — 0.70 and ildm = 0.25, which pro- 
duces a good fit to current CMB data. 

By evaluating the SDSS likelihood we found that 
the couplings are constrained to be 



Q2 < 10^^^ cm^ MeV"^ , 
Qo < 10"^"* cm^ MeV"^ 



(25) 
(26) 



at the 2(7 confidence level in above fiducial cosmology. 
If we compare Eq. (|25|l with Eq. lO, we see that for 
couplings g of order one this bound is saturated if both 
"T-dm and M are of the order of MeV. Smaller values 
of g imply lighter masses for the NIDM and the inter- 
mediate exchanged particle in the scattering process. 
In view of the BBN bound mentioned in Section II, 
nT-dm > 10 MeV, these values are already disfavored, 
so the LSS constraint we obtain is not further con- 
straining the NIDM scenario. In fact, the BBN can 
be weakened if we allows for more exotic features in 
the neutrino density, in particular a neutrino chemical 
potential. Indeed, in this case, the larger contribution 
of DM particles to the Hubble expansion rate resulting 
into a higher ^He mass fraction can be compensated 
by a positive (i.e., more neutrinos than antineutrinos) 
value of the chemical potential;, see, e.g., Ref. p^ . 
Further bounds might be obtained by studying the 
LSS power spectrum at larger wavenumbers, k > 0.2 h 
Mpc~^, taking into account the nonlinear behavior of 
perturbations for very small scales. 

On the other hand, the result for Qq which we recall 
correspond to a intermediate particle and NIDM mass 
degeneracy, is more severely constraining rridm- We 
get in this case rridm > lO/i"*/^ GeV. 



IV. ASTROPHYSICAL CONSTRAINTS 

It is interesting to compare Eqs. (|25|l and (|2(j|l 
with the bounds on Q2 and Qq that can be obtained 
from the propagation of astrophysical neutrinos. The 
most important constraint is provided by observation 
of neutrinos from SN1987A 39], which arc in good 
agreement with the theoretical expectation of neu- 
trino fluxes from type II supernovae. These neutrinos 
have energies of order 10 MeV. The thickness of the 
dark-matter layer that they propagate through is ap- 
proximately / p{l)dl, the integral of the dark-matter 
density along the line of sight I to the LMC. Approx- 
imating the dark-matter density p{l) ~ Po(V^o)^^, 
where po ~ 0.4 GeV cm^"^ is the local density and 
~ 8 kpc our distance from the Galactic center. 
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we find a dark-matter thickness ~ 10^^ MeV cm"^. 
Given the agreement between the predicted and ob- 
served neutrino flux and energy spectrum, we infer 
that neutrinos from SN1987A were not significantly 
absorbed by dark matter along the line of sight, from 
which we get an upper bound ^ 10~^^ cm^ MeV~^ to 
the neutrino-DM opacity for neutrinos of energy ^ 10 
MeV. From this result we obtain the upper bounds 

Q2 < 10 (^[^) cm^Mey-\ (27) 

Qo < 10"^^ cm^ MeV"^ . (28) 

We note that the bound on Q2 is stronger than what 
is obtained using LSS data and of the same order of 
magnitude as the BBN limit corresponding to M > 10 
MeV, while for Qo the stronger bound is still provided 
by Eq. (PI- 

Neutrinos with high energy are likely to be pro- 
duced by a variety of astrophysical sources. Strong 
scattering of these neutrinos off the NIDM when trav- 
eling over cosmological distances of order of tens of 
Mpc implies large energy losses and correspondingly 
a strong deformation of the emitted energy spectrum 
at the source. For light NIDM and intermediate ex- 
changed particles (in the 10-MeV range) the high- 
energy {Ei^ > GeV) scattering cross section behaves 
as adm-u{Eu > lOMeV) g'^/s with s = nidmE^. 
We stress once more that high values for NIDM or 
intermediate F ov U particle mass, though perfectly 
legitimate, implies no observable effects in the LSS 
power spectrum and are thus of no interest for the 
present analysis. Using the definition of Q2 and Qo 
of Section II and the value of the critical density to- 
day, we can evaluate the typical scattering length as 
a function of the neutrino energy as follows 



A. 



1 1 



lOMpC M VLdm Q2 

Xu E„ 1 1 



10 



lOMpc M Qdm Qi 



^'^ c-rn^ MeV^^ , (29) 
lO^^-^cm^ MeV^^ (30) 



If we use the LSS bounds of Eqs. H25|) and H26|l we see 
that for M ~ 10 MeV, the value of \^ is typically very 
large. The effect of interactions with NIDM can only 
affect neutrinos with order-GeV energy over distances 
of 10 100 Mpc assuming the largest value for Q2, 
while the effect is negligible for higher values of Ei, or 
for the Qo case. 

The bounds discussed so far are obtained under the 
assumption that all dark matter is interacting with 
neutrinos. However, if the dark matter is made of 
several components, it is possible that only a frac- 
tion of the dark matter was actually strongly coupled. 
In Fig. ^ we show matter power spectra for a stan- 
dard ACDM model with ^dm = 0.25 and ^tv = 0.7 
and for other models where a fraction of the energy 
density Vlint = 0.7 — JIa of the cosmological constant 
is replaced by interacting dark matter with coupling 
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FIG. 4: Matter power spectra for cold-|-interacting dark 
matter with Q2 = 10"^** cm^ MeV"^ (top panel) and Q2 = 
10"^^cm^ MeV"^ (bottom panel). 



Q2. As we can see, the spectra are quite similar 
if we consider quite large values of Q2, as large as 
10~^*cm^ MeV~^ This is simply due to the fact that 
the interacting component is nearly unclustered on 
large scales, k ^ 0.01ft- Mpc^^. Therefore, adding 
this component or changing the energy density in A is 
nearly equivalent. This degeneracy in the framework 
of NIDM might weaken the current estimates of the 
matter density from galaxy clustering. For higher val- 
ues of Q2 or r^int, oscillations in the power spectrum 
are instead clearly visible. 

In Fig. O we plot constraints on Q2 using SDSS 
P{k) data by allowing this possibility that only a frac- 
tion of the dark matter interacts with neutrinos. The 
overall matter density is fixed at Sim = 0.27 and we as- 
sume a flat universe. As we can see, a smaller ^Inidm 
allows the possibility of relaxing the constraints on Q2. 

Recently, experimental evidence for a "peak" in the 
correlation function of the SDSS luminous red galaxy 
(LRG) distribution at 100 Mpc scales has been re- 
ported 40]. This peak is generally interpreted as the 
imprint of oscillations in the photon-baryon fluid near 
recombination. In the standard CDM framework, this 
peak serves as an absolute ruler and with clustering 
measurements as a function of redshift, one can ex- 
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FIG. 5: Constraints on the Q2 vs Q,nidm plane from SDSS 
P(fc) measurements. An overall matter density of 57m = 
0.27 is assumed with = 0.04. 



low baryon density^. 

As we can see, NIDM is able to mimic the bary- 
onic peak even in the case of a low baryon density 
universe though in the standard CDM, a low baryon 
density will not describe the oscillation. While the 
possibility shown in Fig. is certainly fine tuned, one 
may need to allow for such a scenario when deriving 
constraints on the baryon energy density Q,i,h? from 
large-scale-structure observations if departures from 
the standard description were to be considered. Let 
us note, however, that the above scenario, while mim- 
icking baryonic oscillations in galaxy clustering with 
NIDM, would lead to a different shape for the CMB 
anisotropy power spectrum. A combined analysis is 
therefore a powerful tool for detecting NIDM in mixed 
models and any discrepancy between the value of 
the baryon density derived independently from those 
datasets could hint for an NIDM component. 
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FIG. 6: Degeneracy for the baryon acoustic peak in the 
galaxy correlation function. Three models are plotted: (1) 
a standard CDM model with a baryon density that repro- 
duces the baryon peak (dotted line), (2) a standard CDM 
model with a low baryon density (dashed line), and (3) a 
NIDM model with a low baryon density (solid line). While 
a CDM model with a low baryon density fails to describe 
the peak, a NIDM model with a low baryon density de- 
scribes the peak adequately. The data are taken from the 
LRG analysis of Eisenstein et al. 2005. 



tract strong constraints on the dark-energy time evo- 
lution. However, it is clear that if a portion of the dark 
matter is oscillating as a consequence of strong cou- 
pling with neutrinos, that component may also pro- 
duce an oscillatory behavior and affect the conclusions 
of Ref. • In order to show this effect more quan- 
titatively we compare the LRG data in Fig. with 
the correlation function computed under three model 
descriptions: a standard CDM model with a baryon 
density that reproduces the baryon peak, a combi- 
nation of a standard CDM model with a low baryon 
density, and a model with a NIDM component and a 



V. DISCUSSION AND CONCLUSIONS 

In this paper, we have studied the cosmolog- 
ical consequences of a possible coupling between 
neutrinos and light dark matter with mass in the 
MeV range. We considered two possible behav- 
iors for the thermally-averaged neutrino-DM elastic- 
scattering cross section, either decreasing with tem- 
perature as or constant. We compared the NIDM 
scenario with the large-scale galaxy distribution and 
obtained upper limits on the opacity (ratio of the DM- 
neutrino cross section to the dark-matter mass) of 
O2 < lO-'*^ cm^ MeV-i and Qq < 10-34cm2 MeV'^ 
at the 95% C.L. These limits may be relaxed if one 
consider the possibility that only a fraction of the 
dark matter is made of NIDM. The main cosmologi- 
cal observable for NIDM consists in diffusion-damped 
oscillations in the matter power spectrum. Those 
NIDM oscillations may affect current cosmological 
constraints on neutrinos masses and dark energy from 
galaxy clustering. We have stressed that strongly- 
coupled DM particles would have a non-negligible relic 
abundance today only if an asymmetry between DM 
particle and antiparticle is produced at some early 
stage in the evolution of the Universe, since their den- 
sity would vanish today because of effective annihi- 
lation processes into neutrinos down to temperatures 
much smaller than the DM mass. Detection of NIDM- 
induced oscillations in the LSS power spectrum would 
be a hint for such a non-standard scenario. 



^ As suggested in Ref. |40|. using N-body simulations, to com- 
pare between measurements in the redshift space and predic- 
tions in the real space, we multiply the predicted correlation 
functions by a conversion factor of [1 -1-0.06/(1 -I- (0.06s)^)]^, 
where s is the physical scale measured in Mpc 
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